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Abstract 

We investigate the topological charge in 1+1 dimensional (p^ theory on a lattice with Anti Periodic Bound- 
ary Condition (APBC) in the spatial direction. We propose a simple order parameter for the lattice theory 
with APBC and we demonstrate its effectiveness. Our study suggests that kink condensation is a possible 
mechanism for the order-disorder phase transition in the l-i-l dimensional (^^ theory. With renormalizations 
performed on the lattice with Periodic Boundary Condition (PBC), the topological charge in the renormal- 
ized theory is given as the ratio of the order parameters in the lattices with APBC and PBC. We present a 
comparison of topological charges in the bare and the renormalized theory and demonstrate invariance of 
the charge of the renormahzed theory in the broken symmetry phase. 
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I. INTRODUCTION 



Topology has an important role in nonperturbative quantum field theory on a lattice. For early 
work see, for example, IjJ]. It is of utmost importance in the dual super-conductivity picture of 
quark confinement in QCD In field theories with non-trivial topology, different topological 
sectors are characterized by a conserved topological charge which depends on the boundary be- 
haviour of the theory. This conserved quantity is well-defined in classical continuum field theory. 

In 1-1-1 dimensional classical 0^ field theory, there exists a kink solution with a conserved 
topological charge In the quantum version of the above theory using lattice regularization, 
past studies on kinks focussed on the calculation of the kink mass flQ]. Recently these kinks 
have also been studied in discrete light cone quantization 

In this paper we investigate the topological charge in the context of the ordered (broken sym- 
metry) phase of 1-1-1 dimensional lattice field theory. The short range quantum fluctuations 
lead to the renormalization of the mass and the coupling constant in the topologically trivial sector 
of the theory. On the other hand, long range fluctuations lead to a phase transition in the strong 
coupling region which leads to a breaking of the conservation of topological current and hence the 
vanishing of the topological charge in the disordered (symmetric) phase. We show that the topo- 
logical charge in the renormalized theory remains invariant in the ordered phase. To the best of our 
knowledge, topologocal charge and its invariance in this quantum theory has not been investigated 
before. 

Given that the 0^ theory and the Ising model are in the same universality class, one can also 
define a topolgical charge in the Ising model when considered as a Euclidean quantum field theory 
in 1-1-1 dimensions itzl]. However, our aim in this work is to look at the renormalizations of the 
parameters of the 0^ theory so that a sensible definition of topological charge, invariant in the 
ordered phase of the theory, emerges. 

Condensation of the topological objects has been associated with the mechanism of phase tran- 
sitions in various statistical and quantum field theories |Q. S S (also see works cited in fl). Our 
study indicates that condensation of kink-antikinks is a possible mechanism for the order-disorder 
phase transition in l-i-l dimensional ^'^ theory. 

The plan of this paper is as follows. In Sec. II we present our notation and definition of the 
topological charge in the bare and renormalized lattice theory. Sees. Ill, IV and V are devoted 
to the numerical work with details of various issues that we encounter in the calculation. Finally, 
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Sec. VI contains discussion of our results and conclusion. 



II. TOPOLOGICAL CHARGE AND RENORMALIZATION 



We start from the Lagrangian density in Minkowski space (in usual notation) 



if = i<9^(|)<9>-^mV-^<^>' (2.1) 
which leads to the Lagrangian density in Euclidean space 

Note that in one space and one time dimensions, the field ^ is dimensionless and the quartic 
coupling A has dimension of mass^. 
The Euclidean action is 



5£ = y d^x££E- (2.3) 
Next we put the system on a lattice of spacing a with 

fd^x = a^Y,- (2-4) 

Because of the periodicity of the lattice sites in a toroidal lattice the surface terms will cancel 
among themselves (irrespective of the boundary conditions on fields) enabling us to write 

{d^(l>f = -^dl^ (2.5) 

and on the lattice 

^i^ = ~z [^^+^' + - 2<|)x] . (2.6) 

Introducing dimensionless lattice parameters and Ao by = and Aq = A we arrive at 
the lattice action in two Euclidean dimensions 

^ = - + (2 + f ) I <^>.' + ^ I . (2.7) 

X fl X X 

With antiperiodic boundary condition (APBC), we need to identify an order parameter to char- 
acterize different phases of the theory since (^) is zero in both symmetric and broken phases. In 
previous works on kinks in two dimensional lattice field theory vanishing of the kink mass 



in the symmetric phase was used to identify the critical coupling. Since the calculation of the kink 
mass is rather involved, it is preferable to have a simpler choice. 

For a L? lattice, we propose (^^^^^ = ^[^l-i^^o] order parameter where 

0^ = ( kink g.s kink g.s ) with | kink g.s ) denoting the kink ground state. is computed by 
taking the average of 0x over configurations with APBC in the spatial direction (which effectively 
performs importance sampling around the kink configurations in a Monte-Carlo simulation). 

In the classical theory, the topological charge is given by 



Qclassical = \/ ^diff = J ^diff (2.8) 



B 



where ^diff = 5 [0(°°) —^{—°°)\ and we have used the fact that the mass of the elementary boson 
is MB = \/2m in the broken phase. In the classical theory, the relation between (^diff , mg and A is 
given by 



3ml 

<^>diff=\/^ (2.9) 

which guarantees that the topological charge is +1 (kink sector) and -1 (antikink sector) in the 
broken symmetric phase and zero in the symmetric phase. 

In the lattice theory with APBC in the "spatial" direction, the bare topological charge is defined 

by 



eo = w^<^diff (2-10) 

Due to quantum fluctuations, the classical relation analogous to Eq. ( 12.91 ) between ^^iff, and 
Xq is not obeyed and as a consequence, Qq does not remain +1 or -1 in the ordered phase. 
We propose to define the topological charge in the renormalized theory as 



fi« = J5|?»- ail) 

where = Z being the field renormalization constant. The conventional definition of the 
renormalized coupling Xr in 0^ theory is in terms of the renormalized four-point vertex function. 
The calculation of the four-point vertex function on the lattice is computationally demanding. 
Fortunately, in the broken phase, we can choose a definition lioi of Xr in terms of the renormalized 
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mass iriR and the renormalized vacuum expectation value {^r), determined with periodic boundary 
condition (PBC): 

2 

?iR = 3^. (2.12) 

This definition of involves only the computation of one-point function for (0) and the two point 
function for mr and Z. The details of these computations are provided in Ref. lilliL Using this 
definition of X.r in Eq. (12.1 IL we get 

/-I ^^diff ^diff /O 1 ox 

^^ = 7XT = 7:vr- (2.13) 

The second equality of Eq. (12.131) assumes Z to be the same for ^jiff and ^ and therefore involves 
two unrenormalized quantities, ^^^ff calculated with APBC and (^) calculated with PBC. If 
and (0) are of the same magnitude, the renormalized topological charge is ±1 in the ordered 
phase except for the statistical and systematic errors that occur in the numerical computations of 
the quantities in the numerator and the denominator of Eq. (12.131) . 

ni. CALCULATION WITH PBC : ( ) 

First we discuss the calculation of (0) with PBC. As is well-known, spontaneous symmetry 
breaking cannot occur at finite volume due to the quantum mechanical tunneling phenomena. This 
is due to the fact that at finite volume we have finite degrees of freedom. Tunneling probability 
is, however, inversely proportional to the exponential of the volume and hence is exponentially 
suppressed at large volume lll0lll2l] . Thus in practice, the issue is how large a volume we need for 
the suppression to occur. 

For 0^ theory, tunneling is present for 10^ lattice but is practically absent for 24^ lattice. This is 
demonstrated in Fig. [T] where ^ averaged over all sites for a given configuration is plotted against 
the configuration number. These calculations were done using single hit Metropolis algorithm. 

We find (0) to be sensitive to the various choices of the initial configuration (see Fig. |2l) in the 
simulation using the Metropolis algorithm. For cold starts [^=-1-1 (or -1) for all sites] the selection 
of the vacuum is controlled by the choice of the sign of ^ irrespective of the couplings. On the 
other hand for hot starts (^ = ±1 sprinkled randomly over the entire lattice), the vaccum is picked 
randomly from the two degenerate vacua. 

A more serious problem with the Metropolis algorithm is critical slowing down which makes 
computation near the critical point prohibitively expensive. In the Monte-Carlo simulation, to 
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FIG. 1: Comparison of tunneling in (j)^ theory for FIG. 2: Sensitivity to initial configuration (Metropo- 



couplings mg = —0.5, Ao = 1.8. 



lis algorithm). 



reduce critical slowing down, we have thus resorted to the standard Metropolis update combined 
with a cluster algorithm j3] to update the embedded Ising variables in theory, following the 
method of Brower and Tamayo |14l]. We use Wolff's single cluster variant of the algorithm. This 
method has been used previously in determining the phase diagram of two dimensional lattice 0^ 
theory H . 

In the cluster algorithm however, since the sign of the field of all the members of the cluster 
are flipped in every updation cycle, the algorithm actually enforces tunneling and the configura- 
tion average of ^ is always zero. Thus to get the appropriate nonzero value for the condensate 
we measure (|^|) where ^ = y^j^^^ ^ ^(■^)- To understand the mod let us consider a local order 



sites 



parameter {^{x)). Since the configurations will be selected at random dominantly from the neigh- 
borhood of either vacua in the broken phase, {^{x)) will vanish when averaged over configurations 
thus wiping out the signature of a broken phase. If one uses ( 1 (x) | ) as the order parameter then in 
the broken phase it correctly projects itself onto one of the vacua yielding the appropriate non-zero 
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FIG. 3: (101) versus Ao for various lattice volumes. 
For comparison, classical result is also given. 



FIG. 4: Phase diagram. 



value. The use of this mod, unfortunately, destroys the signal in the symmetric phase completely 
by wiping out the significant fluctuations in sign. However if we choose to use (| yy — E 0(^) !)> 



sites 



it correctly captures the broken phase as well as the symmetric phase. While the sign fluctua- 
tion over configurations are still masked, the fluctuations over sites survive producing ( 1 1 ) = 
correctly in the symmetric phase. 

The volume dependence of ( 1 1 ) is presented in Fig.|3]where we also compare with the classical 
value of the condensate. As can be seen, for small volumes, the signal for phase transition is very 
weak to detect, but 50^ lattice is big enough to observe the transition. 

In Figs. 121 and |3] and all the figures to follow the standard error bars, if not visible, are smaller 
than the symbols unless otherwise stated. 

The phase diagram we obtained for a 5 12^ lattice using PBC is presented in Fig. |4l This agrees 
with the phase diagram obtained in the authors extrapolate their results to infinite 

volume. We observe that our 512^ lattice results are as good as the infinite volume result in Ref. 
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FIG. 5: (j)^ for cold-start(+l) at ml = -0.5 for lat- 
tices (a) 32^, (b) 64^, (c) 128^ and (d) 256^. 

Q. 

Most of our calculations are performed at niQ 
is around 1.95. 



FIG. 6: (j)^ for kink-start at ml = —0.5 for lattices 
(a) 32^, (b) 64^, (c) 128^ and (d) 256^. 

= —0.5. The corresponding critical coupling Xq 



IV. CALCULATION WITH APBC: KINK CONFIGURATIONS AND ^^jg 

Cluster algorithms are known to fail with antiperiodic boundary conditions in theory ll^ . 
For the calculation of ^^iff, because of APBC, we cannot use the cluster algorithm and we resort 
to the standard Metropolis algorithm. In addition to the problems associated with critical slowing 
down, we also face problems associated with computing the profile of an extended object, the 
topological kink in this case and also in fixing its location for the measurement of topological 
charge. 

First we note that the location of the kink cannot be controlled with a cold start as we show 
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in Fig. |5l It is perhaps not easy at first sight, to recognize that most of the configurations seen in 
Fig. 121 are actually kink configurations. To bring it out we refer to Fig. IT] in which we demonstrate 
schematically, how a kink centered near the boundary of a toroidal lattice with APBC on the fields 
should look like. Note that the (L + z) th site is identified with the i th site on a toroidal lattice but 
the sign of the field is flipped on account of APBC on the fields giving it a plateau-like appearance. 
Most of the configurations in Fig. |5lresemble this plateau. 

All such profiles in Fig. |5lare therefore basically kinks or antikinks located near the edge. We 
think that the formation of kinks near the edge are favored by the algorithm because it is clearly 
much easier to generate the plateau-like configurations from a cold start (0 = -fl) compared to a 
kink which involves a flipping of the sign of fields over larger region. Near the critical region of 
course, we expect the formation of kink like configurations to be much easier and it is consistent 
with our observation (Aq = 1 .8 in Fig. |5l). Let us mention at this point that because of translational 
invariance, the position of a kink is of no significance as long as it does not tend to the spatial 
infinities. 

To make our results more transparent we intended to work with kink configurations that are 
centered near the middle of the lattice. The definition of the topological charge that we use (Eq. 
(12.1 11) ) also presupposes that our kink is actually located near the center. 

To obtain such kink configurations we generated configurations by using the kink start (^j^ = — 1 
for < ;c < I and = +1 for | < jc < L) which nicely generates the kink configurations near the 
center (Fig. |6l). 

In passing we remark that the data in both. Fig. |5land Fig. |6lare seen to deteriorate with smaller 
volumes owing to finite size effects which are as usual more serious near the critical point (large 
coupling in this case). We have also observed that there is a small but noticeable increase in the 
size of the kinks (spatial extent over which changes) with an increase of Xq which affects the 
sharpness of the kinks in a small volume near the critical region. 

Reasonably good kink configurations for the measurement of topological charge (with stable 
flat regions) are obtained with the kink start with 128^ and 256^ lattices. 

We have observed kink motion in the vicinity of the critical region. This motion of the kink 
is clearly visible in Fig. [8] where we present for a 512^ lattice for three different time slices. 
Probably because the kink mass is small enough at the couplings, we see this motion. This is 
consistent with the vanishing of the kink mass at the critical point. Incidentally, the equality 
in the magnitude of ^^^^^ among different time slices demonstrates the expected conservation of 
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FIG. 7: The part of the kink between L to 3L/2 gets 
translated to to L/2 (short dash) due to periodicity 
of toroidal lattice and finally gets inverted (solid line 
between to L/2) due to APBC on the field. 



FIG. 8: 0^ for kink start for three different time 
slices showing motion of the kink. 



topological charge. One would have ideally liked to increase the statistics by averaging over the 
time slices at a given x, but the movement of kink makes this not viable in general. However, 
taking into account the fact that the location of the kink is of no consequence in the measurement 
of topological charge, we have actually forced the formation of kinks at the middle (by fixing the 
value of the field at the mid-point in every time slice, to zero) and then taken the average over time 
slices. 



V. KINK CONDENSATION AND TOPOLOGICAL CHARGE 

Let us now come to the discussion of the topological charge. For the use of Eq. (12.131) for the 
topological charge, it is necessary that we get the same phase diagram and critical behaviour for 
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FIG. 9: (101) for PBC and |(^>dig| for APBC versus FIG. 10: Kink condensation shown for 512^ lattice 
Ao at = -0.5. Qo for APBC also plotted. at ml = -0.5. 

the numerator and the denominator of this expression. The simulation results with PBC and APBC 
shown in Fig. I^Jfor 48^ lattice indeed confirm this within our numerical accuracy. Here we would 
also like to remark that Fig. I^Jshows that the disappearance of the kink configuration (characterized 
by I ^jjiff I = 0) and the onset of the phase transition from broken to symmetric phase (characterized 
by (1^1) =0) coincide. It therefore suggests 'kink-condensation' as a possible mechanism for 



the order-disorder phase transition in l-i- 
exist in 2-dimensional Ising model S ^ 



dimensional 0^ theory. Such a connection is known to 
which is believed to be in the same universality class 
and is not totally unexpected in (j)^ theory since it has an embedded Ising variable. Condensation 
of topological excitation has generally been associated with the mechanism of phase transition 
in many statistical and quantum field theories. Kink condensation is clearly visible in Fig. flOl 
where we present for three couplings very close to the critical region. At Xq = 1.95, kink 
configuration is just barely visible, i.e., the boundary value of is still nonzero and passes 
through zero only once. As the coupling increases, the boundary value of ^x reduces to zero 
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FIG. 11: Topological charge calculated with 48^ lat- FIG. 12: Topological charge calculated with 512^ 



tice at ml = —0.5. 
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signalling disappearance of the single kink. However, passes through zero many times showing 
closely packed kink-antikink configurations of varying amplitudes. 

Finally we present the result for the topological charge in Figs. [TT]and[T2l The behavior of the 
topological charge for a range of Ao is depicted in Fig. [Uand Fig. [121 for 48^ and 512^ lattices 
respectively. The figures clearly demonstrate the need for renormalization: It restricts Q^to+l for 
the whole range of Aq investigated in conformity with our expectations. Near the critical region 
the fluctuation of Qr around unity is due to different systematic and statistical errors associated 
with different numerical algorithms used for the calculation of the numerator and the denominator 
of Eq. (12.131 ) especially because the numerator is evaluated with Metropolis algorithm which is 
known to suffer from critical slowing down. We also see by comparing Figs. [TT] and [T^l that there 
is noticeable improvement of data with volume near the critical region. 
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VI. SUMMARY AND CONCLUSIONS 



In this work we have investigated the topological charge in 1+1 dimensional lattice ^ field 
theory. We have shown that with APBC in the spatial direction, lowest energy configuration is a 
kink or an antikink. In order to characterize the different phases of the theory with APBC we have 
proposed a simple order parameter ^^ji^y and we have demonstrated its effectiveness. 

In the process of computing the topological charge we have observed that as the system moves 
on from ordered to disordered phase, the single kink (or antikink) gives way to the occurance of 
a multitude of kink-antikink pairs suggesting kink condensation as a possible mechanism for this 
phase transition. 

A major issue in extracting the topological charge is the renormalization. By a particular choice 
of the renormalized coupling (in the topologically trivial sector), we are able to express the topo- 
logical charge in the renormalized theory as the ratio of renormalized order parameters in the 
lattice theories with APBC and PBC. Provided the field renormalization constants are the same in 
the two cases, the expression for the topological charge in the renormalized theory becomes the 
ratio of unrenormalized order parameters. Making use of this, we have computed the topological 
charge in the renormalized theory and demonstrated that it indeed maintains the value +1 in the 
broken phase. 
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